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Abstract. Consider a weighted and undirected graph, possibly with self- loops, and its corre- 
sponding Laplacian matrix, possibly augmented with additional diagonal elements corresponding to 
the self-loops. The Kron reduction of this graph is again a graph whose Laplacian matrix is obtained 
by the Schur complement of the original Laplacian matrix with respect to a subset of nodes. The 
Kron reduction process is ubiquitous in classic circuit theory and in related disciplines such as electri- 
cal impedance tomography, smart grid monitoring, transient stability assessment in power networks, 
or analysis and simulation of induction motors and power electronics. More general applications of 
Kron reduction occur in sparse matrix algorithms, multi-grid solvers, finite-element analysis, and 
Markov chains. The Schur complement of a Laplacian matrix and related concepts have also been 
studied under different names and as purely theoretic problems in the literature on linear algebra. 
In this paper we propose a general graph-theoretic framework for Kron reduction that leads to novel 
and deep insights both on the mathematical and the physical side. We show the applicability of our 
framework to various practical problem setups arising in engineering applications and computation. 
Furthermore, we provide a comprehensive and detailed graph-theoretic analysis of the Kron reduction 
process encompassing topological, algebraic, spectral, resistive, and sensitivity analyses. Throughout 
our theoretic elaborations we especially emphasize the practical applicability of our results. 
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1. Introduction. Consider an undirected, connected, and weighted graph with 
n nodes and adjacency matrix A G M'^^'^. The corresponding loopy Laplacian matrix 
is the matrix Q G with oflF-diagonal elements Qij = —Aij and diagonal elements 

^ Aij. Consider now a simple algebraic operation, namely the Schur 
complement of the loopy Laplacian matrix Q with respect to a subset of nodes in 
the graph. As it turns out, the resulting lower dimensional matrix Q^ed is again a 
well-defined loopy Laplacian matrix, and a graph can be naturally associated to it. 

This paper investigates this Schur complementation from the viewpoint of alge- 
braic graph theory. In particular we seek answers to the following questions. How 
are the spectrum and the algebraic properties of Q and Qred related? How about the 
corresponding graph topologies and the effective resistances? What is the effect of a 
perturbation in the original graph on the reduced graph, its loopy Laplacian Qred, its 
spectrum, and its effective resistance? Finally, why is this graph reduction process 
and its various properties of practical importance and in which application areas? 
These are some of the questions that motivate this paper. 

Electrical Networks and the Kron Reduction. To show the physical dimen- 
sion and practical importance of the problem setup introduced above, we associate 
an electrical circuit to the graph induced by A. Consider a connected electrical net- 
work with n nodes, branch conductances Aij > 0, and shunt conductances An > 
connecting node i to the ground. By Kirchhoff 's and Ohm's laws the current-balance 
equations / = QV are obtained, where / G W^^^ are the currents injected at the 
nodes, V G W^^^ are the nodal voltages, and the conductance matrix Q G W^^^ 
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is the loopy Laplacian matrix. In various applications of circuit theory and related 
disciplines it is desirable to obtain a lower dimensional electrically-equivalent network 
from the viewpoint of certain boundary nodes (or terminals) a C {l,...,n}, |a| > 2. 
If /5 = {1, . . . , n}\a denotes the set of interior nodes, then, after appropriately labeling 
the nodes, the current-balance equations can be partitioned as 

'Ja 

Gaussian elimination of the interior voltages in equations (1.1) gives an electrically- 
equivalent reduced network with |a| nodes obeying the reduced current-balances 

la + Qac^lB = QredVa , (1-2) 

where the reduced conductance matrix Qred G M'"' ^ I"' is again a loopy Laplacian given 
by the Schur complement of Q with respect to the interior nodes /3, that is, Qred = 
Qaa-Qa^Q^lQpa- The accompanying matrix Qac = -Qa^Q^} ^ mI«|x(^-I«I) maps 
internal currents to boundary currents in the reduced network. In case that //^ is the 
vector of zeros, the (z, j)-element of Qred is the current at boundary node i due to a 
unit potential at boundary node j and a zero potential at all other boundary nodes. 
From here the reduced network can be further analyzed as an |a|-port with current 
injections la + and transfer conductance matrix Qred- 

This reduction of an electrical network via a Schur complement of the associated 
conductance matrix is known as Kron reduction due to the seminal work of Gabriel 
Kron [37], who identified fundamental interconnections among physics, linear algebra, 
and graph theory [33, 38]. The Kron reduction of a simple tree-like network with- 
out current injections or shunt conductances is illustrated in Figure 1.1, an example 
familiar to every engineering student as the Y — A transformation. 



Qo 



Va 



(1.1) 




Fig. 1.1. Kron reduction of a star-like electrical circuit with three boundary nodes 1^, one 
interior node • , and with unit conductances resulting in a reduced triangular reduced circuit. 



Literature Review. The Kron reduction of networks is ubiquitous in circuit 
theory and related applications in order to obtain lower dimensional electrically- 
equivalent circuits. It appears for instance in the behavior, synthesis, and analysis of 
resistive circuits [56, 60, 59], particularly in the context of large-scale integration chips 
[48, 53, 1]. When applied to the impedance matrix of a circuit rather than the admit- 
tance matrix, Kron reduction is also referred to as the "shortage operator" [2, 3, 35]. 
Kron reduction is a standard tool in the power systems community to obtain station- 
ary and dynamically-equivalent reduced models for power flow studies [58, 10, 61], or 
in the reduction of differential-algebraic power network and RLC circuit models to 
lower dimensional purely dynamic models [45, 52, 5, 18, 20]. A recent application of 
Kron reduction is monitoring in smart power grids [17] via synchronized phasor mea- 
surement units. Kron reduction is also crucial for reduced order modeling, analysis. 
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and efficient simulation of induction motors [16, 43] and power electronics [55, 15]. 
Kron reduction is also known in the literature on electrical impedance tomography, 
where Qred is referred to as the "Dirichlet-to-Neumann map" [9, 13, 14]. More gen- 
erally, the Schur complement of a matrix and its associated graph is known in the 
context of (block) Gaussian ehmination of sparse matrices [24, 29, 50, 49], in sparse 
multi-grid solvers [57, 47], and in finite-element analysis [39, 16]. It serves as popular 
application example in linear algebra [25, 23, 54, 40], a similar concept is employed 
in the stochastic complement [41] or the cyclic reduction of Markov chains [8], and a 
related concept is the Perron complement [42, 44] of a matrix and its associated graph. 

This brief literature review shows that Kron reduction is both a practically im- 
portant and theoretically fascinating problem occurring in the reduction of networks 
and their associated matrices. Surprisingly, little is known about the graph-theoretic 
properties of the Kron reduction process. Yet the graph-theoretic analysis of the Kron 
reduction provides novel and deep insights both on the mathematical and the physi- 
cal side of the considered problem. For instance, in the transient stability assessment 
of power networks, it is convenient to analyze the Kron- reduced network and after- 
wards relate its spectral and element-wise connectivity conditions to the non-reduced 
network [18, 20]. In electrical impedance tomography the Kron-reduced matrix Q^ed 
can be directly obtained from boundary measurements or can be constructed from 
effective resistance measurements [51, 22], and the problem is to infer properties of 
Q from Qred [9, 13, 14]. Other graph-theoretic problems occur in the applications of 
smart grid measurements [17] and power flow studies [58, 10], when disturbances such 
as load fluctuations or line outages in the non-reduced network have to be related to 
the Kron-reduced network or vice versa. Also the behavioral analysis of circuits relies 
strongly on Schur complement techniques and methods from graph theory [56, 60, 59]. 
Yet another graph-theoretic problem occurring in numerous applications is minimizing 
the number of links during the reduction process to speed up costly computations as- 
sociated with the network [24, 48, 50, 49, 29, 57, 8, 16, 43, 47, 55, 15, 39]. In summary, 
it might be beneficial to analyze the Kron reduction process from a graph-theoretic 
viewpoint to address the problems above. Finally, the related literature on applied 
matrix analysis [25, 23, 54, 40, 41, 42, 44] shows that the graph-theoretic analysis of 
the Kron reduction process is an interesting mathematical problem in its own right. 

Contributions. The contributions of this paper are two- fold: 
As a first contribution of this paper, we provide a graph-theoretic framework 
and interpretation of the Kron reduction process. This framework encompasses the 
applications of Kron reduction in classic circuit theory [56, 60, 59, 48, 53, 1], electrical 
impedance tomography [9, 51, 22, 13, 14], power flow studies [58, 10, 61], monitoring 
in smart grids [17], and transient stability assessment [45, 52, 5, 18, 20]. We envision 
that our framework for Kron reduction captures various other reduction instances 
of physical problems, it can be applied to more abstract numerical problems as in 
[24, 29, 50, 49, 57, 8, 47], and it can be extended to complex- valued, directed, and 
infinite- dimensional networks, possibly with an underlying dynamical system. 

The second contribution is a detailed analysis of the Kron reduction process. Es- 
sentially, Kron reduction of a connected graph, possibly with self-loops, is a Schur 
complement of corresponding loopy Laplacian matrix with respect to a subset of 
nodes. We relate the topological, the algebraic, and the spectral properties of the re- 
sulting Kron-reduced Laplacian matrix to those of the non-reduced Laplacian matrix. 
Furthermore, we relate the effective resistances in the original graph to the elements 
and effective resistances induced by the Kron-reduced Laplacian matrix. Thereby, we 
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complement and extend various results in the literature on the effective resistance of 
a graph [35, 31, 6, 21]. In our analysis, we carefully analyze the effects of self- loops, 
which model loads and dissipation in the applications of Kron reduction. Finally, 
we present a sensitivity analysis of the algebraic, spectral, and resistive properties of 
the Kron-reduced matrix with respect to perturbations in the non-reduced network 
topology. As a pure mathematical contribution, our analysis of Kron reduction com- 
plements the literature in linear algebra [25, 23, 54, 40] and analogous results on the 
Perron complement side [41, 42, 44]. In our analysis we do not aim at deriving only 
mathematical elegant results but also useful tools for the mentioned applications. We 
believe that our general analysis is a first step towards more detailed results in the 
specific application areas of Kron reduction. 

Paper Organization. The remainder of this section introduces some notation 
recalls some preliminaries in matrix analysis, algebraic graph theory, and concerning 
the effective resistance. Section 2 presents the general framework of Kron reduction 
and reviews various application areas and their specific problem setups. Section 3 
presents the graph-theoretic analysis of the Kron reduction process. Finally, Section 
4 concludes the paper and suggests some future research directions. 

Preliminaries and Notation. Sets: Given a finite set Q, let \Q\ be its cardi- 
nality, and define for n G N the index set = {1, . . . , n}. 

Vectors and matrices: Let Ipxq and Opxg be the p x q dimensional matrices of 
unit and zero entries, and let In be the n-dimensional identity matrix. For vectors, we 
adopt the shorthands Ip = Ipxi and Op = O^xi and define to be vector of zeros of 
appropriate dimension with entry 1 at position i. For a real- valued Id-array {xi}^^^^ 
we let diag({x^}^^^) G M'^^'^ be the associated diagonal matrix. 

Given a real- valued 2d-array {Aij} with i, j G X^, let A G W^^^ denote the 
associated matrix and the transposed matrix. We use the following standard 
notation for submatrices [63]: for two non-empty index sets a^P C let 74[a,y5] 
denote the submatrix of A obtained by the rows indexed by a and the columns 
indexed by /3 and define the shorthands A[a^ f3) = A[a,In \ /3], A{a, = A[In \a,^]^ 
and A{a,/3) = A[Xn \ a, In \ /3]. We adopt the shorthand A[{i}, {j}] = A[iJ] = Aij 
for i, j G Xn, and for x G the notation x[a^ {1}] = x[a] and x(ce, {1}) = x{a). For 
illustration, equation (1.1) can be written unambiguously as 



I[a] 1 




Q[a,a\ 


Q[a, 


a) 


7(a) J 




Q{a,a] 




a) 



V[a] 



If A{a, a) is nonsingular, then the Schur complement of A with respect to A{a^ a) (or 
equivalently the indices a) is the \a\ x \a\ dimensional matrix A/A{a^a) defined by 

A/A{a^ a) = A[a^ a] — A[a^ a)A{a^ a)~^A{a^ a] . 

If A is Hermitian, then we implicitly assume that its eigenvalues are arranged in 
increasing order: Xi{A) < . . . <Xn{A). For a review of matrix analysis we refer to [7]. 

Algebraic graph theory: Consider the undirected, connected, and weighted graph 
G = (Xn^S^A) with node set 1^ and edge set E C x induced by a symmetric, 
nonnegative, and irreducible adjacency matrix A e W^^^ . A non-zero off-diagonal 
element Aij > corresponds to a weighted edge {i^j} G and a non-zero diagonal 
elements An > corresponds to a weighted self-loop {i^i} G £. 

The Laplacian matrix is the symmetric and irreducible matrix defined by L{A) = 
L = diaig{{'^^^-^^ Aij}2=i) — A. Note that self-loops, even though apparent in the 



Kron Reduction of Graphs with Apphcations to Electrical Networks 



5 



adjacency matrix A, do not appear in the Laplacian matrix L. For these reasons and 
motivated by the conductance matrix in circuit theory, we define the symmetric and 
irreducible loopy Laplacian matrix Q{A) = Q = L -\- diagd^d^^} -L^) G R"^^"^. Note 
that adjacency matrix A can be easily recovered from the loopy Laplacian Q as A = 
— Q + diagd^J^-^ - Qij}^^i)^ and thus Q uniquely induces the graph G. We refer to 
Q as strictly loopy Laplacian, respectively as loop-less Laplacian, if the graph induced 
by Q features at least one positively- weighted self-loop, respectively no self-loops. 

For a connected graph ker(L) = span(ln), and all n — 1 remaining non-zero 
eigenvalues of L are strictly positive. Specifically, the second-smallest eigenvalue 
X2{L) is a spectral connectivity measure called the algebraic connectivity. Recall that 
irreducibility of either A, L, or Q is equivalent to connectivity of G, which is again 
equivalent to X2{L) > 0. We refer the reader to [30] for further details. 

Effective resistance: The effective resistance Rij between two nodes i, j G of 
an undirected connected graph induced by a loopy Laplacian Q is defined by 

Rii = {ei - ejfQ\ei - e,) = Q], + Q] • - 2Q\. , (1.3) 

where is the Mo ore- Penrose pseudo inverse of Q. Since is symmetric (follows 
from the singular value decomposition), the matrix of effective resistances R is again 
a symmetric matrix with zero diagonal elements Ru = 0. For two distinct nodes i, j G 
Xn the reciprocal l/Rij is referred to as the effective conductance between i and j. 

REMARK 1.1 (Physical interpretation). If the graph is understood as a resistive 
electrical network with conductance matrix Q, the effective resistance Rij corresponds 
to the potential difference between the nodes i and j when a unit current is injected in i 
and extracted in j. Indeed, in this case the current-balance equations are Ci — ej = QV, 
where V is the vector of potentials with respect to the ground (i.e., nodal voltages). 
The effective resistance Rij defined as the potential difference Rij = (e^ — ej)^V can 
then be obtained via the impedance matrix as Rij = (e^ — ej)^Q\ei — Cj). □ 

The effective resistance Rij can be thought of as a graph-theoretic distance be- 
tween two nodes i and j. In most applications the effective resistance is defined for 
a pure topological network without shunt conductances, that is, a loop-less and uni- 
formly weighted graph with Q = L. We do not restrict ourselves to this case here. We 
refer the reader to [35, 20, 27, 31, 28, 26, 6] for various applications and properties of 
the effective resistance as well as interesting results relating i?, L, Q, ^ and Q~^. 

2. Problem Setup and Applications. 

2.1. The Kron Reduction Process. Consider an undirected, connected, and 
weighted graph G = {In, A) and its associated symmetric and irreducible matrices: 
the adjacency matrix A G M'^^'^, Laplacian matrix L{A), and loopy Laplacian matrix 
Q{A). Furthermore, let a C 2^ be a proper subset of nodes with \a\ > 2. We define 

the (|a| X |a|) dimensional Kron-reduced matrix Qred by 

Qred — Q/Q{a, a) = Q[a, a] - Q[a, a)Q{a, a)~^Q{a, a] . (2.1) 

In the following, we adopt the nomenclature from circuit theory and refer to the nodes 
a and \ a as boundary nodes and interior nodes, respectively. The following lemma 
establishes the existence of the Kron-reduced matrix Qred as well as the structural 
closure property of loopy Laplacian matrices. 

Lemma 2.1 (Structural Properties of Kron Reduction). Let Q G M^><^ 
be a symmetric irreducible loopy Laplacian matrix and let a be a proper subset ofXn 
with \a\ > 2. Then the following statements hold: 
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1) Existence: The Kron-reduced matrix Qred = Q/Q{a,a) is well defined. 

2) Closure: If Q is a symmetric loopy, strictly loopy, or loop-less Laplacian 
matrix, respectively, then Qred 'is a symmetric loopy, strictly loopy, or loop- 
less Laplacian matrix, respectively. 

3) Accompanying matrix; The accompanying matrix Qac — — Q[<^, a)Q{a^ a)~^ 
^ ]^|a|x(n-|Q!|) nonnegative. If the subgraph among the interior nodes is 
connected and each boundary node is adjacent to at least one interior node, 
then Qac 'is positive. If additionally, Q = L is a loop-less Laplacian, then 
Qac = ^ac = —L[a^a)L{a^a)~^ is column stochastic. 

An interesting consequence of Lemma 2.1 is that Qred, as a loopy Laplacian ma- 
trix, induces again an undirected and weighted graph. Hence, Kron reduction, orig- 
inally defined as an algebraic operation in equation (2.1), can be equivalently inter- 
preted as a graph-reduction process^ or as physical reduction of the associated circuit. 
This interplay between algebra, graph theory, and physics is illustrated in Figure 2.1. 




conductance 
matrix Q 

<^ II <^ 

loopy Laplacian 
matrix Q 




KronJJ^reduction 



KronJJ^reduction 



KronJJ^reduction 




transfer conductance 
matrix Qred 

Kron-reduced 
loopy Laplacian Qred 




Fig. 2.1. Illustration of an electrical network with 4 boundary nodes I and 8 interior nodes 
where all line and shunt conductances are of unit value. The associated conductance matrix ( or loopy 
Laplacian) Q and the corresponding graph G are equivalent representations of this network. Kron 
reduction of the interior nodes • results in a reduced electrical network among the boundary nodes 
I with the Kron-reduced matrix Qred ^-^ conductance matrix and the corresponding graph Gred- 



In the following we denote the graph induced by Qred as Gred, and define Ared — 
-Qred + diag({Ej=i,j/iQred[^, and Lred = ^(^red) to be the corresponding 

reduced adjacency and loop-less Laplacian matrices. We remark that Lemma 2.1 is 
partially noted in [5, 25, 56, 60, 54, 57] and present its proof in the following. 

Proof of Lemma 2.1. By definition, Q is (weakly) diagonally dominant since Qu = 
j^i \Qij \ + ^ii foi" all i ^Xn- Due to the irreducibility of Q the strict inequality 
Qii > Xlj=i j^i j^a \Qi3 \ holds at Icast for one i G Xn\a. It follows that Q(a, a) 
is also irreducible, diagonally dominant, and has at least one row with strictly positive 
row sum. Hence, Q(a,a) is invertible [34, Corollary 6.2.27] and statement 1) follows. 

Statement 2) is a consequence of the structural properties of Q and the closure 
properties of the Schur complement [63, Chapter 4], which includes the classes of 
symmetric, positive definite, and M-matrices. Since Q is a symmetric M-matrix, we 
directly conclude that Qred = Q/Q{<^^a) is also a symmetric M-matrix. Hence, Qred 
is a symmetric loopy Laplacian matrix. This fact together with the closure of positive 
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definite matrices under the Schur complement reveals that the class of symmetric 
strictly loopy Laplacian matrices is closed under Kron reduction. To prove the closure 
of symmetric loop-less Laplacians, assume without loss of generality that a = X\a\^ 
and consider the following equality for the row sums of the loop- less Laplacian Q: 



Q[a,a\ 


Q[a,a) 








0|c.| 


Q{a,a] 






In— |q!| 







(2.2) 



Elimination of the second block of equations in (2.2) results in 0|q,| = Qred^\a\^ which 
shows that Qred is a loop-less Laplacian and concludes the proof of statement 2). 

The second block of equations in (2.2) can be rewritten as In-iai = Qac'^\a\i 
which shows that Qac is a column stochastic matrix in the loop-less case. In general, 
Qac = —Q[a^a)Q{a^a)~^ is nonnegative, since —Q[a,a) and the inverse of the M- 
matrix Q{a^ a) are both nonnegative. If additionally each boundary node is connected 
to at least one interior node and the graph among the interior nodes is connected, then 
every row of —Q[a^a) has at least one nonnegative entry and Q{a^a)~^ is positive 
as well (since Q{a,a) is an irreducible non-singular M-matrix [25, Theorem 5.12]) , 
which guarantees positivity of Qac- This completes the proof of statement 3). □ 

The Kron reduction process is ubiquitous in circuit theory and related scientific 
fields. Its general purpose is to construct low dimensional static or dynamic equiva- 
lents of higher dimensional models. In the following we describe different examples of 
Kron reduction in the fields of circuit theory, electrical impedance tomography, power 
flow studies, transient stability assessment, and smart grid monitoring. 

2.2. Kron Reduction in Large-Scale Integration Chips. In circuit theory, 
it is of interest to reduce the complexity of large-scale circuits by replacing them 
with an equivalent circuit with the same terminals (boundary nodes) but with fewer 
branches. This problem occurs in the context of large-scale integration chips [48, 53, 1], 
where the smaller equivalent circuits are used to effectively compute the effective 
resistance between boundary nodes. The circuit reduction problem also stimulated a 
matrix-theoretic and behavioral analysis from the viewpoint of the boundary nodes 
[56, 60, 59, 5] (where boundary nodes correspond to leaves in the modeling framework 
of [60, 59]). For purely resistive circuits the proposed solution in [56, 60, 59, 48, 53, 1] 
to obtain an equivalent circuit is the Kron reduction of the current-balance equations 
/ = QV resulting in the reduced equations (1.2). A particular reduction goal in [48] is 
to reduce the fill-in of the Kron-reduced matrix Qred, see Figure 2.2 for an illustration. 




Fig. 2.2. a) Illustration of an integration chip with a symmetric top-level network connecting 
the interior nodes with the terminals. The Kron reduction of all interior nodes of a network with 
59 terminals results in a Kron-reduced matrix of dimension 59^ with 59^ = 3481 non-zero entries. 
If all but five (specifically chosen) interior nodes are eliminated, then the Kron-reduced matrix is of 
dimension 64^ but has only 1506 non-zero entries. The sparsity patterns of these two Kron-reduced 
matrices are illustrated in subfigures b) and c), which are taken from [48]. 
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The proper choice of the boundary nodes has a tremendous effect on the sparsity of the 
Kron-reduced matrix and saves numerical effort in subsequent computations, which 
is also a pervasive objective in the applications [24, 29, 50, 49, 57, 39, 47, 16]. 

In [48] a constructive algorithm is presented to obtain a sparse Kron-reduced 
matrix with minimal fill-in for computation of the effective resistance. It is argued 
that reduction of a connected and sparse component of Q results in a dense reduced 
component in Qred and the effective resistance among boundary nodes is invariant 
under the Kron reduction process. We remark that these arguments are based on 
numerical observations and physical intuition rather than on mathematical proofs. 
This paper, puts the statements of [48] on solid mathematical ground. We prove 
invariance of the effective resistance under Kron reduction and rigorously show how 
and under which conditions the topology changes from sparse to dense or even com- 
plete. The latter conditions may also be of interest for efficient matrix computations 
[24, 29, 50, 49, 57, 39, 16, 47]. Moreover, our setup encompasses shunt loads and cur- 
rents drawn from the interior network, thereby generalizing results in [56, 60, 59, 5]. 

2.3. Electrical Impedance Tomography. In electrical impedance tomogra- 
phy the goal is to determine the conductivity inside a compact and connected spatial 
domain C from simultaneous measurements of currents and voltages at the 
boundary of Vt^ i.e., from measurement of the Dirichlet-to- Neumann map. Electri- 
cal impedance tomography finds applications in geophysics and medical imaging. A 
natural approach is a discretization of the spatial domain to a resistor network with 
conductance matrix Q [9] . When a unit potential is imposed at boundary node j and 
a zero potential at all other boundary nodes, the current measured at boundary node 
i gives Qred[hj]i ^hc (z, j) element of the reduced conductance matrix (or discretized 
Dirichlet-to-Neumann map) Qred- Other methods iteratively construct the reduced 
impedance matrix (or discretized Neumann-to-Dirichlet map) Q^^^ from measure- 
ments of the effective resistance matrix R [13, 51, 22]. The goal is then to recover the 
original network Q from the reduced network Qred, i-e., inverting the Kron-reduction 
as in Figure 2.3. This is feasible only for highly symmetric networks as considered in 
[13, 14, 9], but generally it is not possible to infer structural properties from Qred to Q. 




Fig. 2.3. In electric impedance tomography the conductivity of the spatial domain 17 is estimated 
by measuring the Kron-reduced matrix Qred boundary ports I. From these measurements the 
conductance matrix Q is re- constructed and serves as a spatial discretization ofQ. 

This paper provides non-iterative methods to relate the effective resistance matrix 
R and the inverse Kron-reduced matrix QI^^ as well as simple explicit formulas to 
relate R and Qred directly for uniform topologies. Furthermore, our analysis allows 
to partially invert the Kron reduction by estimating the spectrum of Q or its effective 
resistance from the spectrum or resistance of Qred- Finally, our framework allows also 
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for dissipation of energy in the spatial domain via discrete loads in the resistor network. 

2.4. Sensitivity of Reduced Power Flow. Large-scale power generation and 
transmission networks are also modeled as interconnected circuits [52, 61]. The nodes 
in a power grid are the 1-connected generators and the buses, see Figure 2.4. Each edge 




Fig. 2.4. Single line diagram, of the New England Power Grid [45], an equivalent schematic 
representation with generators I and buses and the corresponding Kron-reduced network 

connecting two nodes i and j is weighted by a non-zero line admittance Aij = Aji G C 
which is typically inductive, i.e., negative imaginary. Whereas a generator i injects 
current G C into the grid, the load at a bus j draws a constant current Ij G C 
and features a constant shunt admittance Ajj G C. Hence, the power network obeys 
the current-balance equations / = QV, where the nodal admittance matrix Q G C^^"^ 
is the loopy Laplacian induced by the admittances Aij. Alternatively, the current 
balance equations / = QV can be converted to the power flow or power balance 
equations S = F o (5*V*, where o is the entry- wise (Hadamard) product, * denotes 
the conjugate transposed, and S = V o P is the vector of power injections. Depending 
on the application and the set of available input data, it is more convenient to work 
with the linear current-balance equations or the nonlinear power-balance equations. 

A critical task in power network operation is monitoring, control, and optimiza- 
tion of the power flow to guarantee a stable and optimal operating point [61]. Such 
an ideal operating point depends on the phases and magnitudes of the bus voltages 
and the powers and currents injected by the generators. The determining equations 
5* = VoQ* V* are too complicated to admit an analytic solution and often too onerous 
for a computational approach [58, 10]. If a set of boundary nodes a is identified for 
sensing or control purposes and all power injections are set to zero, then all interior 
nodes can be eliminated via Kron reduction leading to the reduced current-balance 
equations (1.2). The corresponding reduced power flow equations are obtained as 
^red = V[a] o Q;,d^[a]*, where ^red = V[a] o /[a]* + V[a] o Q*,/(a)*. 

For both analytic and computational purposes it is important to know the sen- 
sitivity of the reduced power flow w.r.t. changes in the original network topology 
and fluctuations in generation or load. For the lossless case when Q is purely imagi- 
nary, this paper provides insightful and explicit results showing how perturbations in 
weights or topology of Q affect the reduced transfer admittance matrix Qred- We also 
show the effect of shunt loads on the reduced network. In particular, a positive load 
Qii > in the non-reduced network is shown to weaken the mutual transfer admit- 
tances Qredlhj] ill the reduced network and to increase the reduced loads Qred[hi\- 

2.5. Monitoring of DC Power Flow in Smart Grid. The linearized DC 
power flow equations are P = BO^ where P = G is the vector of real power 

injections, G is the vector of voltage phase angles, and B = —^{Q) G is 
the susceptance matrix^ i.e., the loopy Laplacian matrix induced by —^{Aij). The 
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DC power flow often serves as precursor for analytic and computational approaches to 
the nonlinear power flow [61]. Consider now the problem of monitoring an area ^7 of a 
smart power grid with n nodes and equipped with synchronized phasor measurement 
units at the buses a = {ai,a2} bordering the area [17]. Kron reduction of the DC 
power flow P = BO with respect to the interior nodes Tn\a yields the reduced DC 
power flow P[a\ + Bg,cP{(^) = BredO[a], where Pred and Pac are defined analogously 
to Qred and Qac- From here various scalar stress measures over the area Q can be 
defined [17]. Let a G M'"' be the indicator vector for the boundary buses ai, that 
is, (Ji = 1 if z G ai and zero otherwise. The cutset power flow over the area Q 
is Pcut = cr^P[a] + Bg,cP{(^)^ the cutset susceptance is 6cut = cr^^red^^, and the 
corresponding cutset angle is 6>cut = ^cut/^cut = {cr^P[a\ + cr^Pac^(c»^))/(cr^^redcr). 
Hence, the area Q is effectively reduced to two nodes {1,2} exchanging the power flow 
Pcut with angle ^cut over the susceptance bcut^ see Figure 2.5. These scalar quantities 




P[a] P[a]+PacP(c^) ^cut 

Fig. 2.5. Reduction of the area fl with boundary buses M to a single line equivalent {1,2} 
describing the electrical characteristics between the boundary buses ol\ and the boundary buses ol2- 

indicate the stress within the area Vt. For instance, a large cutset angle ^cut could be a 
blackout risk precursor. Of special interest are how load changes, line outages, or loss 
of entire nodes within the area Vt or on its boundary a affect the cutset angle Sent- 

This paper provides a comprehensive and detailed analysis of how changes in topol- 
ogy and weighting of the network affect the Kron-reduced matrix Pred- These results 
include the self-loops in the graph (modeling loads) and can be easily translated to the 
cutset angle ^cut to show its sensitivity with respect to perturbations in the network. 

2.6. Transient Stability Assessment in Power Networks. Transient sta- 
bility in a power network is the ability of the generators to remain in synchronism 
when the network is subjected to large transient disturbances such as faults or loss 
of system components or severe fluctuations in generation or load. The dynamics of 
generator i are typically given by the classic model [52, Chapter 7] 

M-Oi = -A^i + Pm^ - Pe4. ^ ^ , (2.3) 

where {Oi^Oi) are the generator rotor angle and speed, Pm,i > is the mechanical 
power input, and Mi > and P^ > are the inertia and damping constant. The 
active output power injected by generator i into the adjacent bus is Pe,i = 5R(6'i). The 
generator dynamics (2.3) and the power flow equations S = VoQ*^* define the classic 
structure-preserving power network model, an index- 1 differential- algebraic system. 
We refer the reader to [52, 45, 18, 20] for further details and references. 

For the considered static load model, all bus nodes can be eliminated via Kron 
reduction and the network is reduced to the dynamic generator nodes. In the reduced 
network, the generators are coupled via the transfer admittance matrix Qred, and the 
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electrical output power is Pe,i = 5R(^red[^])- For simplicity, we consider the lossless case 
?fi{Q) = Onxn^ which implies 5R(Qred) = 0|q!|x|q!|- Define the effective power input of 
each generator as uji = Pm,i-\-^{Vi X]j=i"' Qld'^^ j]^{a\-\-j)) coupling weights as 

Pij = |^z||^ |^(Qred[^5 j]) > (maximum power transfer between generators i and j), 
then equations (2.3) together with the expression for Pe,i give the swing equations 

MiOi = -DA ^uoi- Y^'^l^ Pij sin(<9i - Oj), ieXn- (2.4) 

In [19], we derived sufficient conditions under which the reduced model (2.4) syn- 
chronizes, i.e., all frequency differences Oi{t) — Oj(t) converge to zero. For notational 
simplicity, we assume uniform damping here, that is, Di = D for all i G a. Then two 
sufficient conditions for synchronization of the model (2.4) are [19] 

\a\mhi{Pij} > rnax {uJi - ujj} , or \2{L{Pij)) > ' {^i-ujjf) ' . (2.5) 

The left-hand sides of conditions (2.5) reflect the connectivity of the weighted graph 
induced by the power transfer Pif the term |a| min^^j{P^j} is a lower bound for 
min^ ^l^'^-^ P^j, the worst coupling of one generator to the network, and X^iLiPij)) 
is the algebraic connectivity of the coupling. The right-hand side of conditions (2.5) 
describes the non-uniformity in effective power inputs uJi in either oo-norm or two- 
norm. In summary, conditions (2.5) can be interpreted as "the reduced network 
connectivity has to dominate the network non- uniformity in effective power inputs." 

For uniformly lower-bounded voltages at all generators > y > (due to local 
control) the analysis of this paper will reveal that the spectral condition on X2{L{Pij)) 
in the reduced network can be converted to the spectral synchronization condition 

HL) > (lZl"Li,,</'^' " ^^•)')'^'^ + max{^red[^,i]} , (2.6) 

where L is the Laplacian of the original lossless power network (weighted by ^{—Aij)) 
and Aped [^5 i\ is the ith load in the reduced network. Similarly, if the effective resistance 
among all generators takes the uniform value R and the effective resistance between 
the generators and the ground is uniform as well, then the results of this paper render 
the element- wise condition on |q^| min^^jjP^j} in the reduced network to a resistive 
synchronization condition in the non-reduced network: 

4 > max {ui - ujj} + max{Ared[i, i]} • (2.7) 

Conditions (2.6) and (2.7) state that the network connectivity has to overcome the 
non-uniformity in effective power inputs and the dissipation by the loads, such that the 
network synchronizes. These conditions bridge the gap from classic transient stability 
analysis [52, 45] to the synchronization analysis in complex networks depending on 
the algebraic connectivity \2{L) [4] or the effective conductance 1/R [36]. 

3. Kron Reduction of Graphs. This section analyzes the algebraic, topologi- 
cal, spectral, and sensitivity properties, as well as the effective resistance of the Kron- 
reduced matrix Qred and its associated graph. Throughout this section we assume 
that Q G M^^^ is a symmetric and irreducible loopy Laplacian matrix (corresponding 
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to an undirected and connected graph with n nodes), and we let a be a proper subset 
of Xn with |a| > 2. For simphcity and without loss of generality, we assume that the 
n nodes are labeled such that a = X\a\- All results in this paper can also be stated 
for an arbitrary labeling of the nodes at the cost of more complicated notation. 

3.1. Preliminary Results on the Augmented Laplacian Matrix and It- 
erative Kron Reduction. The following definitions of the augmented Laplacian 
matrix and the iterative Kron reduction will be central to the subsequent develop- 
ments both for illustration and analysis purposes. 

The role of the self- loops induced by a strictly loopy Laplacian Q G ]R^><"^ can 
be better understood by introducing the grounded node with index n + 1. Then the 
strictly loopy Laplacian Q can be regarded as the principal n x n block embedded in 
the (n + 1) X (n + 1) dimensional augmented Laplacian matrix 



Q 


-diag{{Au}2^,)ln ■ 




V" A-- 



(3.1) 



where A G ]R^><^ is the adjacency matrix corresponding to Q. The augmented Lapla- 
cian Q is the Laplacian of the augmented graph G with node set V = {X^, n + 1} and 
edge set S = ^augment}- Here a node z G 2"^ is connected to the grounded node 
n + 1 via a weighted edge {i, n + 1} G ^augment if and only if An > 0. A graph G and 
its associated augmented graph G are illustrated in Figure 3.1. 



G ■ o oo oo G 

oo oo oo 
^ CO o no 




Fig. 3.1. Illustration of the graph G associated with the electrical network from Figure 2.1 and 
the corresponding augmented graph G with additional grounded node (). 

This augmentation and embedding process is inspired by Fiedler's work [26] on non- 
singular M- matrices. The following lemma presents some interesting properties of Q. 

Lemma 3.1 (Properties of augmented Laplacian). Consider the symmetric 
and irreducible strictly loopy Laplacian Q G W^^^ and the corresponding augmented 
Laplacian matrix Q G R("'+i)x("'+i) defined in (3.1). The following statements hold: 

1) Algebraic properties: Q is an irreducible and symmetric loop-less Lapla- 
cian matrix. ^ 

2) Spectral properties: The eigenvalues of Q and Q interlace each other: 

= Ai(Q) < Ai(Q) < A2(Q) < A2(Q) < • • • < Xn{Q) < K{Q) < K+i{Q) • 

3) Kron reduction: Consider the strictly loopy Laplacian Qred = Q/Q{<^^ and 
the loop-less Laplacian Q red — Q/Q{{c^i'f^-^^}i{(^i'f^-^^}) 7 both obtained by Kron 
reduction of the interior nodes Xn\oi. Then the following diagram commutes: 
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^ auerment ^ 

Q — Q 



Kron reduction 
of Xn \ a 



Kron reduction 

of Xn \ a 



^^^^ augment ^ Qred 
In equivalent words, Qred 'is the augmented Laplacian associated to Qred; ^-^v 



Q 



red 



Qred 


-diag({Ared[i,i]}iea)l 


-1^' dmg{{Ared[i:i]}iea) 





(3.2) 



Properties 2) and 3) of Lemma 3.1 intuitively illustrate the effect of self-loops on 
the spectrum of Q and its Kron-reduced matrix. Specifically, the elegant relationship 
3) implies that the reduced self-loops A^edi^,^] can be equivalently obtained by the 
last column or row of Qred- In short, Kron reduction can equivalently be applied to 
the strictly loopy network G or to the augmented loop-less network G. 

Proof of Lemma 3.1. Property 1) follows trivially from the construction of the 
augmented Laplacian Q. Property 2) is a direct application of the interlacing theorem 
for bordered matrices [34, Theorem 4.3.8], where = Ai(Q) < Ai(Q) since Q is an 
irreducible loop-less Laplacian and Q is non-singular. In property 3), the (1,1) block 
of the matrix on the right-hand side of identity (3.2) follows by writing out the Schur 
complement of a matrix partitioned in 3 x 3 blocks, as in the proof of the Quotient 
Formula [63, Theorem 1.4]. The other blocks follow immediately since Kron reduction 
of the loop- less Laplacian Qred yields again a loop- less Laplacian by Lemma 2.1. □ 

Gaussian elimination of interior voltages from the current- balance equations / = 
QV can either be performed via Kron reduction in a single step, as in equation (1.2), 
or in multiple steps, each interior node i G {l,...,n — \a\} at a time. The following 
concept of iterative Kron reduction addresses exactly this point. 

Definition 3.2 (Iterative Kron reduction). Iterative Kron reduction as- 
sociates to a symmetric irreducible loopy Laplacian matrix Q G R^^^ and indices 
{1, . . . , \a\}, a sequence of matrices Q^ G R(^-^^^('^-^\ ^ G {1, . . . , n— |a|}, defined by 

Q' = Q'-'/QU, (3.3) 

where Q^ = Q and ki = n-^l—£, that is, Qi~^^ is the lowest diagonal entry of Q^~^ . 

If the sequence (3.3) is well-defined, then each Q^ is a loopy Laplacian matrix 
inducing a graph by Lemma 2.1. Before going further into the details of iterative 
Kron reduction, we illustrate the unweighted graph corresponding to Q^ (the sparsity 
pattern of the corresponding adjacency matrix) in Figure 3.2. The following observa- 
tions can be made: (1) The connectivity is maintained. (2) At the ^th reduction step 
a new edge between two nodes appears if and only if both were connected to ki before 
the reduction, and (3) all other edges persist. (4) Likewise, a new self- loop appears 
at a node i k£ ii and only if i was connected to ki and ki featured a self-loop before 
the reduction. The next subsection turns these observations into rigorous theorems. 

In components, Q^ is defined according to the celebrated Kron reduction formula 
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1) 



CHD— «: 






Fig. 3.2. Sparsity pattern (or unweighted graph/ topolgical evolution) corresponding to the iter- 
ative Kron reduction (3.3) of a simple graph with 3 boundary nodes I and 7 interior nodes • . The 
dashed red lines indicate the newly added edges in a reduction step. Steps 4) - 6) are omitted. 



which nicely ihustrates the step-wise Gaussian ehmination. For a weh-defined se- 
quence we let and be the corresponding adjacency and loop-less 
Laplacian matrix of the ^th reduction step. The following lemma states some impor- 
tant properties of iterative Kron reduction. In particular, the iterative Kron reduction 
is well-posed and ultimately results in the Kron-reduced matrix. 

Lemma 3.3 (Properties of iterative Kron reduction). Consider the ma- 
trix sequence defined via iterative Kron reduction in equation (3.3). The 
following statements hold: 

1) Well-posedness: Each matrix , £ e {1, . . . ,n — \a\}, is well defined, and 
the classes of loopy, strictly loopy, and loop-less Laplacian matrices are closed 
throughout the iterative Kron reduction. 

2) Quotient property: The Kron-reduced matrix Qred = Q/Q{(^^Q^) can be 
obtained by iterative reduction of all interior nodes e In \ ol, that is, 
Qred = Q"'"'"'. Equivalently, the following diagram commutes: 

Kron reduction of \ 



Qred — Q 



-\a\ 



Q'^Q^^ . gn- 

iterative Kron reduction 



-|a| + l 



J 



3) Diagonal dominance: The ith row sum of , i G {1, . . . ^n — i}, is given by 



Qii ~ ^ii 




a: 



iki y 



1 



1 

kike 



= 0, 
>0. 



(3.5) 



Property 3) of Lemma 3.3 confirms the closure of strictly loopy and loop-less 
Laplacians and further reveals that the weights of the self-loops are non-decreasing. 

Proof of Lemma 3.3. Statement 2) is simply the Quotient Formula [63, Theorem 
1.4] stating that Schur complements (or Gaussian elimination for that mater) can 
be taken iteratively or in a single step. Furthermore, the Quotient Formula states 
that all intermediate Schur complements exist for^G{l,...,n — |a|}. This fact 
together with the closure properties in Lemma 2.1 proves statement 1). For notational 
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simplicity and without loss of generality, we prove statement 3) for ^ = 1 and ki = n. 
Note that = A, = = and consider the ith row sum of given by 

n— 1 n— 1 / 

j=i j=i V 

n-1 

where we used equality (3.4), the identities Q = L-\- didig{{Aii}'^^^) , ^^=i =^ii + 

Ain^ and Yl^=i ^jn = ^nn- Siucc Ann ^ (duc to property 1) nonnegative row sums 
follow also in the general case), we are left with the two cases stated in (3.5). □ 

3.2. Topological, Spectral, and Algebraic Properties of Kron Reduc- 
tion. In this section we finally begin our characterization of the properties of Kron 
reduction. We start by discussing how the graph topology of G, or equivalently the 
sparsity pattern of Q, changes after Kron reduction of the interior nodes. 

Theorem 3.4 (Topological Properties of Kron Reduction). Let G, Gred, 
and G he the undirected weighted graphs associated to Q, Qred = Q/Q{<^:<^), (^'^d the 
augmented loopy Laplacian Q, respectively. The following statements hold: 

1) Edges: Two nodes i^jea are connected by an edge in G^ed ^/ ctnd only if 
there is a path from i to j in G whose nodes all belong to {i^j} U {In \ a). 

2) Self- loops: A node i E a features a self-loop in G^ed ^/ ^^c? only if there 
is a path from i to the grounded node n + 1 in G whose nodes all belong to 
{i} U {Xn \ a). Equivalently, a node i E a features a self-loop in G^ed ci'^d 
only if i features a self-loop in G or there is a path from i to a loopy interior 
node j eln\cy whose nodes all belong to {i^j} U (2"^ \ a). 

3) Reduction of connected components: // the interior nodes (3 C Xn\a 
form a connected subgraph of G, then the boundary nodes a C a adjacent to 
P in G form a clique in Cred- Moreover, if one node in {3 features a self-loop 
in G, then all boundary nodes adjacent to {3 in G feature self-loops in Gred- 

The topological evolution of the graph corresponding to the iterative Kron reduc- 
tion (3.4) - the sparsity pattern of A^ - is illustrated in Figure 3.2. Statement 1) of 
Theorem 3.4 can be observed in each reduction step, statement 2) is visible in the 
third step, and statement 3) is visible in the final step of the reduction in Figure 3.2. 

In the loop-less case. Theorem 3.4 is also proved in [5, Theorem 4.20 and Theorem 
4.23] and in [25, Theorem 14.2]. In the strictly loopy case. Theorem 3.4 can be ob- 
tained from [40, Lemma 3.1]. Given our prior results on iterative Kron reduction and 
the augmented Laplacian matrix, the proof of Theorem 3.4 is rather straightforward. 

Proof of Theorem 3.4- To prove Theorem 3.4 we initially focus on the loop- less 
case and the reduction of a single interior node k via the one-step iterative Kron 
reduction (3.4). Due to the closure of loop-less Laplacian matrices under iterative 
Kron reduction, c.f. Lemma 3.3, we can restrict the discussion to the non-positive 
off-diagonal elements of = Q/Quk inducing the mutual edges in the graph. Any 
non-zero and thus strictly negative element Qij is rendered to a strictly negative 
element Q]j since the first term on the right-hand side of equation (3.4) is strictly 
negative and the second term is non-positive. Therefore, all edges in the graph induced 
by Qij persist in the graph induced by Q\- . According to the iterative Kron reduction 
formula (3.4), a zero element Qij = is converted into a strictly negative element 



Qin Qj 
Qnn 



n-1 



L^nn ~l~ ^nn 



L^nn ~l~ ^nn 



^ ^ i^jn) — ^ii H" ^in 



J = l 



^nn ^ -^nn 



An 
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Qjj < if and only if both nodes i and j are adjacent to k. Consequently, a reduction 
of node k leads to a complete graph among all nodes that were adjacent to k. In 
summary, statement 1) holds in the loop- less case and for a single reduction step. 

Recall from Lemma 3.3 that the one-step reduction of all interior nodes is equiva- 
lent to iterative reduction of each interior node. Hence, the arguments of the previous 
paragraph can be applied iteratively, which proves statement 1) in the loop- less case. 

The results in statement 2), the strictly loopy case, follow simply by applying 
the same arguments to the augmented Laplacian Q defined in (3.1). Alternatively, an 
element-wise analysis of Aj^ together with equation (3.5) leads to the same conclusion. 

Finally, statement 3) follows directly by applying statement 1) and 2) to the 
connected component (3. This completes the proof of Theorem 3.4. □ 

By Theorem 3.4, the topological connectivity among the boundary nodes be- 
comes only denser under Kron reduction. Hence, the algebraic connectivity \2{L) - 
as a spectral connectivity measure - should increase accordingly. Indeed, for the graph 
in Figure 3.2 (with initially unit weights), we have \2{L) = 0.30 < A2(i^red) = 0.45. 
Physical intuition also suggests that loads in an electrical circuit weaken the influ- 
ence of nodes on another. Thus, self-loops should weaken the algebraic connectivity 
A2(^red) the rcduccd network. In the following we confirm these intuitions and 
present some facts on the spectrum of the Kron-reduced matrix. 

Theorem 3.5 (Spectral Properties of Kron Reduction). The following 
statements hold for the spectrum of the Kron-reduced matrix Qred = Q /Q{oL^a): 

1) Spectral interlacing: For any r G X\a\ it holds that 

Xr{Q) < KiQred) < A.(Q[c.,a]) < A,+,_|,|(0) . (3.6) 

2) Effect of self- loops: Specifically, for any r G T\a\ it holds that 

Ar(^red) + max{Ared[^, i]} > K{L) + min{AiJ , (3.7) 
K{Lred) + min {Ared[i,i]} < Xr+ri-\a\{L) + max{Aii} . (3.8) 

In the loop-less case, inequality (3.6) implies non-decreasing algebraic connectiv- 
ity \2{L) < A2(i^red)- Foi" instance, the graph in Figure 2.1 with zero- valued self-loops 
satisfies \2{L) = 0.39 < A2(i^red) = 0.69. However, in the strictly loopy case inequali- 
ties (3. 7)- (3. 7) imply that self-loops weaken the algebraic connectivity tremendously. 
The same graph with unit- valued self-loops satisfies X2{L) = 0.39 > A2 (I/red) = 0.29. 

To prove Theorem 3.5, we recall the interlacing property for Schur complements. 

Lemma 3.6. (Spectral Interlacing Property, [23, Theorem 3.1]). Let A be 
a Hermitian positive semidefinite matrix of order n and let /3 he a, non-empty proper 
subset of In such that is nonsingular. Then for any r G 2'^-|/3|, it holds that 

Xr{A) < Xr{A/A[p,p]) < Xr{A{p,p)) < A,+|^|{A). 

Proof of Theorem 3.5. Since Q is a loopy Laplacian matrix and hence positive 
semidefinite. Lemma 3.6 can be applied with j3 = Xn\oL and results in the bound 
(3.6) in statement 1). To prove statement 2), we recall Weyl's inequality [34, Theorem 
4.3.1] for the spectrum of the sum of two Hermitian matrices A, B of dimension n: 

\k{A) + Ai(5) < \k{A + 5) < \u{A) + \n{B) , keXn. (3.9) 

In order to continue, let r G X|q,| (or equivalently, r G a) and consider the following 
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set of spectral equalities and inequalities resulting in the spectral bound (3.7): 

KiLred) = KiQred - diag({ Aed ^] > K{Qred) - UiaX^^c. { Aed ^] } 

> XriQ) -m^^i^a{Aed[hi]} = A^(L + diag({AiJ^^i)) - maXiea{^red[^, ^]} 

> Xr{L) -\-mmi^Xr,{Aii} - maxi^a{Aed[i,i]} , 

where we subsequently made use of the identity Lred = Qred — diagd^dredi^, ^ll^ea)? 
Weyl's inequality (3.9), the fact Ai(- diag({Ared[^, ^]}iea)) = - maxiea{^red[^, ^]}, 
Lemma 3.6, the identity Q = L -\- diag({Aii} -L^), and again Weyl's inequality (3.9) 
with Xi {disig{{ Aii}2^-^^)) = mmi^Xni^u}- The bound (3.8) follows analogously. □ 

The following theorem summarizes some algebraic properties of the Kron reduc- 
tion process. In particular, it quantifies the topological properties stated in Theorem 
3.4, it quantifies the reduced self- loops occurring in Theorem 3.5, and it states that 
the edge and self-loop weights among the boundary nodes are non-decreasing, as seen 
in Figure 2.1. Furthermore, Theorem 3.7 confirms the intuition that the class of undi- 
rected connected graphs is closed under Kron reduction due to closure of irreducibility, 
and it reveals some more subtle properties concerning the effect of self-loops. 

Theorem 3.7 (Algebraic Properties of Kron Reduction). Consider the 
Kron-reduced matrix Qred Ci'^d the accompanying matrices Qac = —Q[o:, o:)Q{o:, o:)~^ 
and Lac = —L[a^ a)L{a^ a)~^ defined in Lemma 2.1. The following statements hold: 

1 ) Closure of irreducibility: // Q is irreducible, then Qred '^s irreducible. 

2) Monotonicity of elements: For all i^j G a it holds that Qred[hj] ^ 
(5[a,a][z,j]. Equivalently, A^edl^j] ^ j] for alli^j G a. 

3) Effect of self- loops I: Define Ai = An > 0, for i ^X^, so that loopy and 
loop-less Laplacians Q and L are related by Q = L -\- diag({A^}^^x^)- Then 
the Kron-reduced matrix takes the form 

Qred = L/L{a, a) + diag({AJ,ea) + S , (3.10) 

where S = L^c{In-\a\ + diag({AJ,^2:^\^)L(a, a)-^)-Miag({AJi^x^\c.)^rc 
is a symmetric nonnegative \a\ x |a| matrix. Furthermore, the reduced self- 
loops satisfy Ared[i^ z] = A^ + Ej=l''' Qeic[iJ]^\a\+j foriea. 

4 ) Effect of self-loops II: // the subgraph among the interior nodes Tn\a is 
connected, each boundary node a is connected to at least one interior node 
in Xn\ a, and at least one of the interior nodes has a positively weighted 
self-loop, then S and Qac ^^^e both positive matrices. 

Statements 1) and 2) are not surprising given our knowledge from Theorems 3.4 
and 3.5. Statement 3) reveals an interesting fact that can be nicely illustrated by con- 
sidering the reduction of a single interior node k with a self- loop A/^ > 0. In this case, 
the matrix S in identity (3.10) specializes to the symmetric and nonnegative matrix 
S = Ck- L{k, k]L[k, k) e M(--i)x(n-i)^ ^i^ere = Ak/{Lkk{Lkk + A^)) > 0. Hence, 
the reduction of node k decreases the mutual coupling {i^j} in Q/Qkk by the amount 
Ck ' Aik Ajk > and increases each self- loop i in Q/Qkk by the corresponding amount 
Ck • Aik Aik > 0. This argument can also be applied iteratively. The complete reduc- 
tion of all nodes in statement 4) implies, under quite general topological conditions, 
that a single positive self-loop in the interior network will affect the entire reduced 
network by decreasing all mutual connections and increasing all self- loops weights. 

For the proof of Theorem 3.7, we recall the following identities for the inverse of 
the sum of two matrices: 
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Lemma 3.8. (Sherman-Morrison Formula, [32]). Let A,B e C'''' and let A 
be non-singular, li A-\- B is nonsingular, then 

{A^ B)-^ = A-^ - A-\I ^ BA-^)-^BA-K (3.11) 

If additionally B = A- uv^ for A G M and u,v eW, then 

(A + A • uv^)-^ = A-^ - -^--^A-VA-^ . (3.12) 

Proof of Theorem 3.7. In the loop-less case, the spectral inequality (3.6) in Theo- 
rem 3.5 implies non- decreasing algebraic connectivity A2 (I/red) ^ ^2{L) > and thus 
irreducibility of Lred- In the strictly loopy case, note that the Kron- reduced graph 
features the same edges (excluding self- loops) as in the loop- less case, by Theorem 
3.4. Thus, connectivity and irreducibility follow, which proves statement 1). 

The element-wise bound Qredlhj] ^ ceJi^^i] follows directly from [12, Lemma 
1], where this bound is stated for the reduction of one node. By Lemma 3.3, a one-step 
reduction is equivalent to iterative one- dimensional reductions. Hence, [12, Lemma 
1] can be applied iteratively and yields Qredlhj] < Q[<^5 ce][^5 j]- For the off-diagonal 
elements i 7^ j, this bound is readily converted to 74red[^7 j] ^ ^[<^7 Qf] [^7 j]- The same 
bound follows for the diagonal elements since diagonal dominance is non-decreasing 
under Kron reduction, see Lemma 3.3. This completes the proof of statement 2). 

Identity (3.10) in statement 3) follows by expanding the Kron- reduced matrix as 

Qred = Q/Q{a,a) = (L + diag({AJ^exJ)/(iv(a,a) +diag({AJ^exAa)) 

= diag({A Jiea) + ^[q^, c^] - L[<^, cx) {L{a, a) -h diag({AJ^ex^\c,))"^L(a, a] 

= L/L{a, a) + diag({AJ,ea) + 5' , 

where S is defined statement 3), and we applied the matrix identity (3.11) in the fourth 
inequality with A = L(a, a) and B = diag({A^}^^x„\Q;)- By Lemma 3.3, the one-step 
Schur complement Q/Q{a^a) is equivalent to iterative one-dimensional reduction of 
all interior nodes I^i \ the matrix = Q^~^ /Q^'^ki reduction step 

is again a loopy Laplacian. If we abbreviate the self-loops at the ^th reduction step 
by Af = Af-, then can be reformulated according to identity (3.10) as 

= Q'-^/Qi-l = L'-'/Li-l + diag({A|}--/) + , (3.13) 

where = C£ - L^~^{k£^k£]L^~^[ki^k() = q • A^~^(ki^ki]A^~^[ki^ki) is a symmetric 
and nonnegative matrix and q = A^^/ {Lk^ki{Lk^ke_ + ^kg)) — ^- Iterative application 
of this argument implies that S must also be a symmetric and nonnegative matrix. 

To obtain an explicit expression for the reduced self- loops, re-consider the identity 
(2.2) defining the self-loops of L. In the general loopy case identity (2.2) reads as 
Qln = A. Block- Gaussian elimination of the interior nodes yields Qredl|a| = A[a] -h 

(5ac^(ce)- Hence, the ith row sum of Qred satisfies 74red[^5^] = Yl^j^iQred[h j] = 

+ Z]j=l"' Qac[hj]^\a\-\-j- This Completes the proof of statement 3). 
Under the assumptions of statement 4) , the positivity of Qac follows from Lemma 
2.1. To prove positivity of 5*, note that iterative reduction of all but one interior 
node yields one remaining interior node = h. According to equality (3.13), 

reduction of this last loopy node gives the matrix = Ch ■ A^{h^ h]A^[h^ h). Under 
the assumptions of statement 4), Theorem 3.4 implies that h features a self-loop and is 
connected to all boundary nodes. It follows that c/^ > and > for all i G 
Therefore, Sh is a positive matrix, and the same can be concluded for S. □ 
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3.3. Kron Reduction and Effective Resistance. The physical intuition be- 
hind the Kron reduction and the effective resistance in Remark 1.1 suggests that the 
transfer conductances Qred[hj] related to the corresponding effective conductances 
1 / Rij . The following theorem gives the exact relation between the Kron-reduced ma- 
trix Qred^ the effective resistance matrix i?, as well as the augmented Laplacian Q. 

Theorem 3.9 (Resistive Properties of Kron Reduction). Consider the 
Kron-reduced matrix Qred = Q/Q{a^a), the effective resistance matrix R, and the 
augmented Laplacian matrix Q defined in (3.1). The following statements hold: 

1 ) Invariance of under Kron reduction: The effective resistance Rij between 
any two boundary nodes i^jeais equal when computed from Q or Qred-' 



Rij 



2) Invariance under augmentation: If Q is a strictly loopy Laplacian, then 
the effective resistance Rij between any two nodes i^j G In is equal when 
computed from Q or Q: 

Rij = {ei - ej)^Q~^{ei - ej) = {ci - ej)^Q\ei - ej) , ij G In • (3.15) 

In other words, statements 1) and 2) imply that, if Q is a strictly loopy Laplacian, 
then the following diagram commutes: 

Q 



augment 



Kron reduction 
oiXn\a 



Qred 




Kron reduction 
of \ a 



Qred 



augment 

3) Effect of self-loops: If Q is a strictly loopy Laplacian and Ri 



{ei 



ej)^L\ei — Cj), i^j ^In, is the effective resistance computed from the corre- 
sponding loop-less Laplacian L, then for all i^j G In it holds that Rij < Rij. 
Figure 3.3 illustrates the ideas of Theorem 3.9 with six simple graphs. 



Go 



6z 



G, 



O^ll ^33 O 



Gl,red [B- 



A.d[l,3] 



Go 




G. 



^23 



Fig. 3.3. Illustration of Theorem 3.9 with six graphs: According to statement 1), the effective 
resistance R13 between the boundary nodes H is equal when computed in the graph Gi or in the 
corresponding Kron-reduced graph Glared- According to statement 2), the effective resistance Ris is 
equal when computed in the strictly loopy graph G2 (resp. G3) or in the corresponding augmented 
loop-less graph G2 (resp. G3) with additional grounded node (). According to statement 3), the 
effective resistance R13 in the strictly loopy graphs G2 and G3 is not larger than in the loop-less 
graph Gi (with equality for the pair {Gi, G2} and strict inequality for the pair {Gi, G3} ). 

Identity (3.14) states that the effective resistances between the boundary nodes 
are invariant under Kron reduction of the interior nodes. Spoken in terms of circuit 
theory, the effective resistance between the terminals a can be obtained from either 
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the impedance matrix or the transfer impedance matrix QI^^- Identity (3.15) gives 
a resistive interpretation of the self-loops: the effective resistance among the nodes 
in a strictly loopy graph G is equivalent to the effective resistance among the corre- 
sponding nodes in the augmented loop-less graph G with additional grounded node. 
The commutative diagram in Theorem 3.9 combines the invariance identity (3.14) 
with the self-loop identity (3.15). According to statement 3), the self-loops do not 
increase the effective resistance. This is in accordance with the physical interpretation 
in Remark 1.1 when viewing the self-loops as shunt conductances. 

To prove Theorem 3.9, we need the following identities relating R and L. 

Lemma 3.10 (Alternative Derivation of the Effective Resistance). Let 
^ ^ j^nxn ^ symmetric irreducible loop-less Laplacian matrix. Then for any S ^ 

(L + (S/n) Inxn)"' = L^^ (l/Sn) l^xn . (3.16) 

Consider the matrix R of effective resistances defined in (1.3). For 6 it holds that 

Rij = {ci-ejY L\ei-ej) = {ci-ej)^ {L-\-{S/n) lnxn)~^{ei-ej) , ij eXn- (3.17) 

If further n > 3, then by taking node n as reference the following identity holds: 

Rij = {ei - ej)^L\ei - ej) = (e^ - Cj)^ L{n,n)~^{ei - ej) , ij G Xn-i • (3.18) 

Proof Identity (3.16) is a simple generalization of [31, Theorem 5]. Since 
Inxnlnxn = n • l^xn and LL^ = L = In - (1/n) • Inxn (by definition of 
via the singular value decomposition), identity (3.16) can be verified since the prod- 
uct of the left-hand and the right-hand side of (3.16) equal the identity matrix. The 
identity (3.17) follows then by multiplying equation (3.16) from the left by (e^ — ej)^ 
and from the right by (e^ — Cj). Let L := L(n, n), then it follows from [27, Appendix 
B, eq. (17)] that L~^^ = L\.-L\^-L]^-L\^. The identity (3.18) can now be directly 
verified. For an alternative derivation of identity (3.18), we refer to [28]. □ 

Proof of Theorem 3.9. We begin by proving statement 1) in the strictly loopy case 
when Q is invertible and = . Note that we are interested in the effective resis- 
tances only among the nodes a, i.e., the |a| x |a| block of . The Schur complement 
formula [63, Theorem 1.2] gives the |a| x \a I block of as (Q/Q(a,a))-i = Q"^. 
Consequently, for i^j^a the defining equation (1.3) for the effective resistance Rij 
is simply rendered to Rij = (e^ — ^j)^ Q^^i^i — ^j), the claimed identity (3.14). 

In the loop-less case when Q = L is singular, a similar line of arguments holds 
on the image of L. Let ^ > and consider the modified and non-singular Laplacian 
L = L -\- {S/n)lnxn' Due to identity (3.16) we have that = -\- {l/6n)lnxn- 
Since {ci — ej)^lnxn{^i — ^j) = 0, we can rewrite identity (3.17) in expanded form as 

Rij = i^i - ejf{L^ + {1/Sn)lnxn){ei - Cj) = (e^ - ejfL-^a - Cj) . (3.19) 

As before, the \a\ x |a| block of L is given by L-^ as {L/L{a,a))-^ [63, Theorem 
1.2]. Consequently, the identity (3.19) is for i^jea rendered to 

Rij = {ei - ejf{L/L{a,a))-\ei - Cj) . (3.20) 

Note that the right-hand side of (3.19), or equivalently (3.20), is independent of S 
since the matrices are evaluated on the subspace orthogonal to 1^, the nullspace 
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of L as (5 I 0. Thus, on the image of L the hmit of the right-hand side of (3.20) 
exists as (5 ^ 0. By definition, L^^ acts as regular inverse on the image of L, and 
equation (3.20) is rendered to Rij = (e^ — ej)^{L/L{a^a))'^{ei — ej). FinaUy, recaU 
that Lred = L/L{a^a) which yields the claimed identity (3.14) in the loop-less case. 

To prove statement 2), note that the strictly loopy Laplacian Q is invertible (due 
to irreducible diagonal dominance [34, Corollary 6.2.27]). Hence, the defining equation 
(1.3) for the resistance features a regular inverse. The matrix Q can also be seen as the 
principal nxn block of the augmented Laplacian Q, that is, Q = Q(n + 1, n + 1). The 
identity (3.15) follows then directly from identity (3.18) (with n replaced by n + 1). 

To prove statement 3), we appeal to Rayleigh's celebrated monotonicity law and 
short/cut principle [21]. Since the Laplacian L induces the same graph as Q with node 
n + 1 cut off, the monotonicity law states that the effective resistance Rij in the graph 
induced by L is not smaller than the effective resistance Rij in the graph induced by 
Q. The latter again equals the effective resistance in the graph induced by Q due to 
identity (3.15). Equivalently, for i, j G 2^ it holds that Rij = (e^ — Cj)^ L\ei — Cj) > 
(ci - CjYQ^Ci - Cj) = (ci - ej)^Q~^(ei - Cj) = Rij. □ 

Theorem 3.9 allows to compute the effective resistance matrix R from the transfer 
impedance matrix Q^ed* interested in a converse result. Iterative methods 

constructing QI^^ from R can be found in [13, 51, 22]. It is also possible to recover the 
(pseudo) inverse of the loopy Laplacian Q, the augmented Laplacian Q, or the corre- 
sponding Kron-reduced Laplacians directly from the effective resistance matrix R. 

Lemma 3.11 (Impedance and Effective Resistance Matrix). Let Q G M^^^ 
be a symmetric irreducible loopy Laplacian matrix. Consider the following three cases: 

1) Loop-less case: Let R G R^^^ be the matrix of effective resistances. Then 
for i^j G In the following identity holds: 

1 ^ -j^ n n 

Q^[i,j] = L^[i,j] = ~-(Rij --Y,{Rik + Rik) + -^J2J2^ke) ■ (3.21) 

^ \ TL Til / 

k=l k=li=l 

2) Strictly loopy case: Consider the grounded node with index n -\- 1, the 
corresponding augmented Laplacian matrix Q G r(^+^)><(^+i) defined in (3.1)^ 
and the corresponding matrix of effective resistances R G M(^+^)>^(^+^) defined 
in (1.3). Then the following two identities hold: 

^ ^ n+l ^ n+1 n+1 

k=i ^ ^ k=i i=i 

ij ein+i, (3.22) 

Q~' [h j] = \ {R[h n + 1] + R[j, n + 1] - i?[i, j]) , i, j G In (3.23) 

3) Kron reduced case: The identities (3.21)^ (3.22)^ and (3.23) also hold 
when , Q\ and on the left-hand sides are replaced by Q\^^, Q\q^, and 
Qred^ 'respectively, and n on the right-hand sides is replaced by \a\. 

Proof. Identity (3.21) is stated in [6, Theorem 4.8] for the weighted case and in 
[62, Theorem 7] for the unweighted case. According to statement 2) of Theorem 3.9, 
the resistance is invariant under augmentation. Hence, identity (3.21) applied to the 
augmented Laplacian Q yields identity (3.22). Identity (3.23) follows directly from 
[6, Theorem 4.9] or [26, Theorem 2.1]. According to statement 1) of Theorem 3.9, the 
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effective resistance is invariant under Kron reduction, and thus the effective resistance 
corresponding to Qred is simply R[a^ a]. Hence, the formulas (3.21), (3.22), and (3.23) 
can be applied to the Kron-reduced matrix Qred as stated in 3). □ 

The formulas stated Theorem 3.9 and Lemma 3.11 allow to compute the effective 
resistance matrix R in the original non-reduced network from the (pseudo) inverse 
of the Kron-reduced Laplacian Qred, and vice versa. In some applications, such as 
the transient stability problem presented in Subsection 2.6, it is desirable to know 
an explicit algebraic relationship between R and Qred without the (pseudo) inverse. 
However, such an explicit relationship between R and Qred can be found only when 
closed- form solutions of Q\q^^ Q^ed ^Ld known. These are generally not avail- 
able. Generally, it also seems infeasible to relate bounds on R to bounds on Qred since 
element-wise bounding of inverses of interval matrices is known to be NP-hard [11]. 

Fortunately, closed forms of the inverses of QI^^^ Q^ed ^Ld derived in 

an ideal electric network, where the effective resistances among the boundary nodes 
are uniform, as well as the effective resistances between the boundary nodes and the 
ground. In this case, the following theorem states the equivalence of such an ideal 
network and uniform transfer conductances in the Kron-reduced network. 

Theorem 3.12 (Equivalence of Uniformity in Effective Resistance and 
Kron Reduction). Consider the Kron-reduced Laplacian Q^ed = Q/Q{ol^ a) and the 
corresponding adjacency matrix A^-^d- Consider the following two cases: 

Loop- less case: Let R G M^^^"^ he the matrix of effective resistances defined in 
(1.3). Then the following two statements are equivalent: 

1) The effective resistances among the boundary nodes a are uniform, i.e., there 
is r > such that Rij = r for all distinct i^j G a; 

2) The weighting of the edges in the Kron-reduced network is uniform, i.e., there 
is a > such that Ared[hj] = a > for all distinct i^j G a. 

If both statements 1) and 2) are true, then it holds that r = 

Strictly loopy case: Consider additionally the grounded node n + 1 and the 
augmented Laplacian matrices Q and Q^ed defined in (3.1) and (3.2)^ respectively. 
Let R G M(^+i)><(^+i) be the matrix of effective resistances in the augmented network. 
Then the following two statements are equivalent: 

3) The effective resistances among the boundary nodes a are uniform, i.e., there 
is r > such that Rij = r for all distinct i^j G a. The effective resistances 
between all boundary nodes a and the grounded node n-\-l are uniform , i.e., 
there is g > such that R[i^ n -\- 1] = g for all i G a; 

4) The weighting of the edges and the self-loops in the Kron-reduced network is 
uniform, i.e., there are a > and b > such that Aredlhj] = cl > and 

^red 

\i,i] = b>0 for all distinct ij G a. 
If both statements 3) and 4) cire true, then it holds that r = |^|^^^ and g = 

REMARK 3.13 (Physical networks and graph topologies with uniform resistance). 
The assumption of uniform effective resistances in statements 1) and 3) is by no 
means artificial but rather corresponds to an ideal network, where all boundary are 
electrically uniformly distributed in the network, with respect to each other and with 
respect to the shunt loads. In the applications of electrical impedance tomography 
and smart grid monitoring, this assumption can be met by choosing the boundary 
nodes corresponding to sensor locations. In the transient stability problem in power 
networks the generators corresponding to the boundary nodes are distributed over the 
power grid ideally in such a way that the loads can be effectively balanced. Otherwise 
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the power grid is vulnerable to voltage collapse and frequency instabilities. Hence, 
the uniformity assumptions are ideally met in man-made networks. 

Independently of engineered electrical networks, the following examples show that 
uniform resistances among a set of boundary nodes a occur for various graph topolo- 
gies. Here we neglect self- loops and edge weights introducing additional degrees of 
freedom to construct graphs with uniform resistances. In the trivial case, \a\ = 2, 
Theorem 3.12 reduces to [35, Corollary 4.41] and the effective resistance among the 
a nodes is clearly uniform. Second, if the boundary nodes are 1-connected leaves of a 
highly symmetric graph among the interior nodes, such as a star-shaped tree, a com- 
plete graph, or a combination of these two, then the effective resistance among the 
boundary nodes is uniform. Third, the effective resistance in large-scale small-world 
networks is uniform among sufficiently distant nodes [36]. Fourth, with increasing 
number of nodes the effective resistance in random geometric graphs converges to a 
degree-dependent limit [46], which is uniform for various geometries and node dis- 
tributions. Fifth and finally, geometric graphs such as lattices and their fuzzes are 
special random geometric graphs. According to the previous arguments, the resistance 
among sufficiently distant lattice nodes becomes uniform in the large limit. □ 

In order to prove Theorem 3.12, we need the following Laplacian identities. 

Lemma 3.14 (Uniform Laplacian Matrices and their Inverses). Let a > 
and b >0 and consider the loopy (n x n)- dimensional Laplacian matrix Q = a{nln — 
Inxn) + ^In Corresponding to a complete graph with n nodes, uniform positive edge 
weights a > between any two distinct nodes, and nonnegative and uniform self-loops 
b>0 attached to every node. The following statements hold: 

1 ) For zero self-loops b = 0, the pseudo inverse of Q is the loop-less Laplacian 



1 



1 



Q — 2 ' (j^-^'Ti Inxn) • 



2) For positive self-loops b> 0, the inverse of Q is the positive matrix 



Q- 



- {nin 



b{an + b) 

3) Consider the augmented Laplacian Q defined in (3.1) and given as 

Q 



a{nLn - In) + bin 


-bin ' 




n ■ b 



The pseudo inverse of Q is given by the (augmented) loop-less Laplacian 



c(nln - In) + din 


-din ' 


[ -di'n 


n • d 



where d = l/(6(n-h 1)^) and c= -d- {a - (n -h 2)6)/(an -h 6). 

Proof Lemma 3.14 can be verified by direct computation. Statement 2) follows 
since Q - = - Q = In- Statements 1) and 3) follow since Q and (resp. 
Q and Q^) satisfy the Penrose equations [7, Theorem 6.1.1]. Hence, the loop- less 
(augmented) Laplacian (resp. Q^) is the unique Moore- Penrose inverse. □ 

We now have all three ingredients to prove Theorem 3.12: the invariance formulas 
(3.14)-(3.15) for the effective resistance stated in Theorem 3.9, the relations between 
effective resistance and the Kron-reduced impedance matrix in statement 3) of Lemma 
3.11, and the Laplacian identities in Lemma 3.14. Given these formulas, the proof of 
Theorem 3.15 reduces to mere computation. For the sake of brevity, it will be omitted. 
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3.4. Sensitivity of Kron Reduction to Perturbations. In the final section 
of our analysis of Kron reduction we discuss the sensitivity of the Kron-reduced matrix 
Qred to perturbations in the original matrix Q. A number of interesting perturbations 
can be modeled by adding symmetric matrix W eW^^'^ and considering the perturbed 
loopy Laplacian matrix Q = Q + W ^ where Q is the nominal loopy Laplacian. 

The case when 1^ is a diagonal matrix is fully discussed in Theorem 3.7. A 
perturbation of the form when is a non-zero matrix and all other entries of 

W are zero can model the emergence, loss, or change of a self- loop or an edge among 
boundary nodes. In this case, the perturbation acts additively on Q^ed = Q/Q{<^^ as 

Qred = Q/Q{a. a) = Qred + W[a, a] . (3.24) 

If the perturbation affects only the interior nodes, then W{a, a) is a non-zero matrix 
and all other elements of W are zero. Inspired by [5, 61], we put more structure on 
the perturbation matrix W and consider symmetric rank one perturbations of the 
form W = A ■ {ci — ej){ei — e^)-^, where A G M. Such a perturbation changes the 
weight of the edge {i, j} in the network from Aij to Aij -h A and can model the loss 
or the appearance of a new link {i, j} with weight A^j = A. Since a perturbation 
among the boundary nodes is fully captured by (3.24), the more interesting case of a 
perturbation among the interior nodes is presented in the following theorem. 

Theorem 3.15 (Perturbation of the Interior Network) . Consider the Kron- 
reduced matrix Q^q^ = Q/Q{a,a), the accompanying matrixQ^c = —Q[cy^ a)Q{a, a)~^ 
and a symmetric rank one perturbation^ W = A • (e^+|Q,| — ej+|c^|)(e^+|Q,| — ej_^\a\)^ 
for distinct i^j G X^-i^i such that the perturbed matrix Q = Q -\- W remains an 
irreducible loopy Laplacian. The following statements hold: 

1 ) Algebraic perturbation: Q^ed undergoes the rank one perturbation 

1 + A • Rint[l,j\ 

where Rint[i^j] — (e^ — ej)^Q{a^ a)~^{ei — Cj) > 0. 

2) Resistive perturbation: Let R and R be the matrices of effective resis- 
tances corresponding to Q and Q, respectively. For any k,i ein holds that 

Rm = Rm - , . ||(e/c - ei)^Q\ei^oc\ - ej+\cx\)\\l • (3.26) 

// A > (respectively A <{)) then Rki < Rm (respectively Rm > Rm)- 
Note that the term Rint[hj] statement 1) can be interpreted as the effective 
resistance between the perturbed nodes in the interior network induced by (5(a,a). 
Likewise, the physical interpretation of the term Qg,c{ei — ej)A{ei — ej)^ = QacW{a^ a) 
is well-known in network theory [61]. The perturbation W has the same effect on the 
network equations / = {Q-\-W)V as the current injection / = —WV^ that is, the per- 
turbation of the interior edge {i, j} by a value A is equivalent to injecting the current 
A • {Vi_^\a\ ~ ^■+|q!|) iiito the jth interior node and extracting it from the ith interior 
node. In the reduced network equations (1.2) the current injection / translates to the 
current injection Qac-^(ce) = —Qaci^i — ~ ^i)'^^(^) i^^o the boundary nodes. 

Elimination of V{a) then yields I[a] = (5red^[<^]- Finally, the additive term in iden- 
tity (3.26) resembles the sensitivity factor in network theory [17, 61]. From Remark 



^The perturbation W affects the edge between the ith and the jth interior node by a value A. 
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1.1, we see that {ek — e£)^Q^ {ei_^\(^\ — e^+ic^l) corresponds to the potential drop between 
nodes k and ^ if a unit current is injected in the ith interior node and extracted at 
the jth interior node. As before, the current flowing along the perturbed edge {i, j} 
in the interior network is redistributed in the perturbed network in identity (3.26). 

Starting from identity (3.25) various spectral bounds can be derived for the rank 
one perturbations of Q and Qred- For instance, for A < Weyl's inequalities (3.9) give 



Ar(Qred) > K{Qred) > K{Qred) + . 5 ' ||Qac(ei - ej)\\2 , V G X|q,| . 

These bounds can be further related to Q and Q via the interlacing inequalities (3.6) 
or [34, Theorem 4.3.4]. Finally, we remark that the proof of Theorem 3.15, presented 
in the following, can also be applied to more complex symmetric perturbations. 

Proof of Theorem 3.15. Since the perturbed matrix Q = Q+VK is a symmetric and 
irreducible loopy Laplacian, the reduced matrix Qred = Q/Q{a^a) exists by Lemma 
2.1. By employing the matrix identity (3.12), the Schur complement Qred simplifies to 

Qred={Q + W)/{Q{a,a) + W{a,a)) = {Q ^W) / {Q{a,a) ^ A{ei - ejf {ei - ej)) 
Q[a,a)Q{a,a)~^{ei - ej){ei - ej)^Q{a,a)~'^Q{a,a] 



= Qred + A 



This equality further simplifies to identity (3.25) with R^^t as defined above. This 
concludes the proof of statement 1). For the proof of statement 2), we initially consider 
the strictly loopy case. In this case, can be obtained from identity (3.12) as 

1 + A(e^+|c,| - ej^\a\VQ~^{ei^\a\ - ^j+\a\) 

Multiplication of from the left by {e^ — ^i)^ and from the right by {e^ — e^) yields 
identity (3.26). In the loop-less case when Q is singular, the same arguments can be 
applied on the image of Q by considering the non-singular matrix Q + ((^/n) l^xn for 
^ 7^ and identity (3.17). This results in the more general identity (3.26). The second 
part of statement 2) follows then again from Rayleigh's monotonicity law [21]. □ 

4. Conclusions. We studied a graph-theoretic analysis of the Kron reduction 
process. Our analysis is motivated by various example applications spanning from 
classic circuit theory over electrical impedance tomography to power network appli- 
cations. Prompted by these applications, we presented a detailed and comprehensive 
graph-theoretic analysis of the Kron reduction process. In particular, we carried out 
a thorough topological, algebraic, spectral, resistive, and sensitivity analysis of the 
Kron-reduced matrix. This analysis led to novel results in algebraic graph theory and 
new physical insights in the application domains of Kron reduction. We believe our 
results can be directly employed in the application areas of Kron reduction. 

Of course, the results contained in this paper can and need to be further refined 
to meet the specific problems and graph topologies in each particular application 
area. Our analysis also demands answers to further general questions, such as the 
extension of this work to directed graphs [25] or complex-valued Laplacian matrices 
occurring in all disciplines of circuit theory [1, 53, 58, 52]. The results on the effective 
resistance should be extended to the related but more general network sensitivity 
factors [17, 61]. The presented topological analysis should be rendered to a synthesis 
of optimal reduction patterns in terms of sparsity similar to [47, 48, 24, 29, 49, 50]. 



26 



F. Dorfler and F. Bullo 



Some of our results have also analogous results for the Perron complement [42, 44] 
occurring in Markov chain applications [41], but many results still have to be mapped 
from Perron complements to Schur complements and vice versa. Finally, it would be 
of interest to analyze the effects of Kron reduction on centrality measures, clustering 
coefficients, and more general graph-theoretic metrics than the effective resistance. 
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